CeFePO: A Heavy Fermion Metal with Ferromagnetic Correlations 
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The ground state properties of CeFePO, a homologue of the new high temperature superconduc- 
tors RFePnOi- x F x , were studied by means of susceptibility, specific heat, resistivity, and NMR 
measurements on polycrystals. All the results demonstrate that this compound is a magnetically 
non-ordered heavy fermion metal with a Kondo temperature Tk ~ 10 K, a Sommerfeld coefficient 
7 = 700 mJ/molK 2 and a mass enhancement factor of the order of 50. Analysis of the susceptibility 
data and of the spin relaxation time indicate that the strong electronic correlation effects originate 
from the Ce-4/ electrons rather than from Fe-3d electrons. An enhanced Sommerfeld- Wilson ra- 
tio R = 5.5 as well as a Korringa product So/TiTK 2 ~ 0.065 well below 1 indicate the presence 
of ferromagnetic correlations. Therefore, CeFePO appears to be on the non-magnetic side of a 
ferromagnetic instability. 

PACS numbers: 75.20.Hr, 76.60.-k, 71.20.Eh 



The compounds RTPnO (R: rare earth, T: transition 
metal, Pn: P or As) have started to attract considerable 
attention because of the recent discovery of superconduc- 
tivity with a transition temperature T c exceeding 50 K 
in the flFeAsOi-^ series [J 0, H Q, being the high- 
est T c except for cuprate systems. This high T c super- 
conducting state appears to be of unconventional nature 
and related to the disappearance of a Fe-based spin den- 
sity wave. While recent reports focus on the importance 
of electronic correlation effects due to 3d-electrons close 
to a magnetic state, the homologous compounds with R 
= Ce are attractive candidates for 4/-electron induced 
strong correlation effects. Thus last year, we presented a 
detailed study of the properties of CeRuPO and CeOsPO 
[E S 0] and demonstrated that the former one is a rare 
example for a ferromagnetic (FM) Kondo lattice with a 
FM Ce ordering temperature Tc —15 K and a Kondo 
temperature Tk ~ 10 K, while the latter one shows anti- 
ferromagnetic (AFM) order of stable trivalent Ce below 
TV = 4.5 K. More recently, superconductivity at 41 K 
was also reported in CeFeAsOi-^Fz [J], but the relevance 
of 4/-correlations could not yet be analyzed because of 
the very limited amount of experimental data. In this 
report we present a complete study of the basic physical 
properties of CeFePO using susceptibility x(T), specific 
heat C(T), and resistivity p(T) measurements as well as 
NMR as a local probe. The analysis of our results indi- 
cates that CeFePO is a magnetically non-ordered heavy 
fermion metal with strong 4/-based electronic correla- 
tions and FM fluctuations, while Fe correlations seem to 
be of minor importance. After submission of this pa- 
per, a very precise investigation of LaFePO was reported 
by T.M. McQueen et. al. @. These authors conclude 
that LaFePO is a non-magnetic metal with only weak 
exchange enhancement from spin fluctuations. 

Polycrystalline samples were synthesized using a Sn- 
flux method in evacuated quartz tubes as described in 
[H, 0] . Several powder X-ray diffraction patterns recorded 
on a Stoe diffractometer in transmission mode using 
monochromated Cu-K a radiation (A = 1.5406 A) con- 



firmed the P4-/nmm structure type and the formation 
of single phase CeFePO. The lattice parameters a — 
3.919(3) A and c = 8.330(5) A refined by simple least 
squares fitting were found to be in good agreement with 
the reported structure data jlOj . Susceptibility measure- 
ments were performed in a commercial Quantum Design 
(QD) magnetic property measurement system MPMS. 
p(T) was determined down to 0.4 K using a standard AC 
four-probe geometry in a QD physical property measure- 
ment system PPMS. The PPMS was also used to mea- 
sure the specific heat C(T) with a standard heat-pulse 
relaxation technique. NMR measurements are carried 
out by monitoring the spin-echo intensity as a function 
of applied field at a fixed frequency (12.1 MHz, 27 MHz, 
76.4 MHz, and 130 MHz). Shift values are calculated 
from the resonance field H* by K = (H L - H*)/H* 
whereas the Larmor field Hr, is given by using H3PO4 
with 31 K = as a reference compound. Relaxation rates 
31 (l/Ti)(T) are determined by using the saturation re- 
covery method. 

We start by briefly presenting the results of x{T) of 
CeFePO. Below 350 K, xi T ) of CeFePO increases with 
decreasing temperature following a Curie- Weiss law down 
to 100 K (Fig. [3] inset (b)). The effective moment /x e ff = 
2.56 fiB, being almost identical to the expected value 
/i c ff = 2.54 hb for a free trivalent Ce ion, and the mod- 
erate Weiss temperature Oyy — -52 K indicate a nearly 
stable trivalent Ce state and the absence of a signifi- 
cant contribution of a paramagnetic Fe-moment in this 
temperature range. Below 100 K, x(T) deviates from 
the Curie- Weiss law, increasing stronger with decreas- 
ing temperature than expected, which is attributed to 
crystal electric field (CEF) effects. Below 10 K, the 
measured bulk susceptibility becomes strongly field de- 
pendent, showing a pronounced increase for low fields 
fj,oH < 1 T, while for larger fields x(T) levels off at a 
constant value xo = 660 x 10 _9 m 3 /mol. While such a 
field dependence is usually attributed to paramagnetic 
defects or impurities, the NMR results suggest this to 
be an intrinsic property of CeFePO (see below). No 
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FIG. 1: (Color online) 4/-increment to the specific heat of 
CeFePO plotted as C 4f /T vs T in different magnetic fields. 
Inset (a) shows p{T) and (b) p(T 2 ). 

anomaly corresponding to magnetic order could be re- 
solved at any field. Thus, a preliminary analysis of these 
susceptibility data suggests a paramagnetic Kondo lat- 
tice. This conclusion is confirmed by the results of the 
specific heat measurements, plotted as C 4 ? /T vs T in 
Fig. [TJ The 4/-contribution C 4 ? to the specific heat 
was obtained by subtracting the reported specific heat 
data of LaFePO @|. Below 10 K, C 4f /T increases log- 
arithmically down to 1 K, where it levels off at a con- 
stant value 7 = 700 mJ/molK 2 . No anomaly indicating 
a phase transition can be resolved above 0.4 K. Appli- 
cation of a small magnetic field pqH < 1 T does not 
change C 4 f{T), while a larger field leads to a smooth 
decrease of C 4 f /T at low temperatures. The former re- 
sult confirms the absence of long-range magnetic order, 
while the latter one is incompatible with spin-glass-type 
behavior. Therefore, the large C 4 ' /T value at low T has 
to be attributed to heavy fermions and 7 corresponds to 
a strongly enhanced Sommcrfeld coefficient in the elec- 
tronic (heavy- fermion-based) specific heat. Comparing 
this with the Sommerfeld coefficient of the nonmagnetic 
LaFePO (7 = 12.5 mJ/molK 2 ) Q, one estimates a low- 
temperature enhancement factor of the order of 50 for 
the heavy quasiparticles in CeFePO. The magnetic en- 
tropy S{T) was calculated by integrating the measured 
C 4f (T)/T over T for T < 10 K. The entropy S(T) 
reaches 0.5i?ln2 at 5 K, indicating a Kondo temperature 
Tk ~ 10 K. The absolute values of C(T), its temperature 
dependence as well as the deduced Kondo temperature 
are very close to those of the archetypical heavy fermion 
metal CeCu2Si2 [ll|. The temperature dependence of the 
electrical resistivity p{T) shown in the inset (a) of Fig. [1] 
presents the typical behavior of a Kondo lattice system 
and thus confirms the above conclusion. Below 300 K, 
p(T) decreases almost linearly, before dropping rapidly 
below 30 K due to the onset of coherent Kondo scatter- 
ing. Below T = 1 K (see inset (b)), p(T) follows a p oc T 2 
dependence indicating the formation of a Landau-Fermi 
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FIG. 2: (Color online) 31 P field-sweep NMR spectra at 
v — 76.4 MHz and different temperatures. Vertical dashed 
lines indicate the Larmor field obtained from a reference com- 
pound. Inset shows a typical powder spectrum, together with 
the simulation at T = 50 K (arrows indicate resonance fields 
H* for the H _L c and H || c direction). 



liquid ground state in agreement with the specific-heat 
data. However, the absolute value of p(T) is rather high; 
very likely, this is not intrinsic of CeFePO but due to the 
granularity of the polycrystalline sample. 

In order to get a deeper insight into the many-body 
phenomena in this compound, we used 31 P-NMR as a mi- 
croscopic local probe. A typical set of spectra at different 
temperatures is shown in Fig. [2] for a NMR frequency of 
76.4 MHz (corresponding to 4.47 T). One single narrow 
31 P-NMR line as expected from the crystal structure was 
found at room temperature (right panel of Fig. [5]) . It de- 
velops strong anisotropy towards lower temperatures (left 
panel of Fig. The shape is characteristic of a powder 
pattern from a spin 1=1/2 nucleus in a tetragonal sym- 
metry. The broadened 31 P-NMR lines could be simulated 
consistently at all temperatures with shift-tensor compo- 
nents K a i)(T) and K C (T) corresponding to the H _L c 
and H \\ c directions, respectively (inset in Fig. [5J). The 
relative shift values 31 K a b(T) and 31 K C (T) are both posi- 
tive, in contrast to the results obtained for CeRuPO and 
CeOsPO |U. While 31 K C (T) is almost independent of 
T, 31 K a b(T) increases strongly with decreasing tempera- 
tures (Fig. 0]). Above 15 K, 31 K iso = (2K ab /3) + (K c /3) 
tracks the bulk susceptibility x(T) obtained by MPMS 
measurements at 5 T (inset (a)). Below 10 K, the line 
shape of the spectra becomes field dependent. While at 
high fields (poH > 6 T) the line shape observed at higher 
temperatures is maintained down to 2 K, reducing the 
field below 4 T leads to a strong broadening of the contri- 
bution related to H _L c and to a shift of its maximum to 
lower K values. In contrast, the contribution connected 
with H | c is not affected. Since foreign phases or impu- 
rities usually do not contribute to the NMR signal 
and since impurities would not show such an anisotropic 
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FIG. 3: (Color online) Shift component K a t as a func- 
tion of temperature in different magnetic fields as indicated 
(Field values were calculated from NMR frequency using 
( 31 7nmr/27r) = 17.10 MHz/T). Inset (a) shows the isotropic 
shift 31 iiTi S o = (2K a b/3) + (K c /3) versus isotropic (bulk) sus- 
ceptibility from MPMS measurements (5 T). Inset (b) shows 
X _1 (r), solid line corresponds to a Curie-Weiss (CW) fit. 



behavior, this broadening cannot be explained by disor- 
der or impurities. Such a broadening is typical for the 
onset of short-range ferromagnetic correlations. These 
correlations are likely to be responsible for the field de- 
pendent low T increase in x(T), too. Since only the com- 
ponent H _L c is affected in the NMR spectra, only the 
basal-plane components of the magnetic moments start 
to be ferromagnetic correlated. This strong anisotropy 
supports that these correlations originate from the Ce- 
moments. A fit of the whole spectra still lead to a rea- 
sonably well defined K a b value, but the resulting K is 
now field dependent. At high fields, K a i, increases down 
to the lowest T, while at low fields it shows a clear maxi- 
mum. This specific field dependence appears to indicate 
that CeFePO is close to a ferromagnetic instability. 

Finally, we present spin-lattice-relaxation (Ti) data on 
CeFePO. T\ measurements were performed as a function 
of temperature at two different frequencies, 27 MHz and 
76.4 MHz (corresponding to 1.59 T and 4.47 T, respec- 
tively), at the H _L c position of the anisotropic NMR 
line (left arrow in inset of Fig. [3] The nuclear magnetiza- 
tion recovery curves M(t) could be fitted at any temper- 
ature and field with a standard single exponential func- 
tion. 31 (l/TiT) plotted in Fig. 0] increases monotonously 
over two orders of magnitude upon cooling down from 
300 K to 4 K. At lower temperatures, 31 (TiT)~ 1 stays 
constant at about 20 (sK) _1 . Since NMR data on the 
reference compound LaFePO are presently not available, 
we note that the homologue LaRuPO has a 31 (TiT)~ 1 
value ~ 1 (sK) _1 , more than one order of magnitude 
smaller than that of CeFePO. Even 1/TjT of 75 As in 
doped LaFeAOi-^ [H 03 is smaller than that of 
31 P in CeFePO indicating that the relaxation is domi- 



FIG. 4: (Color online) (1/TiT) as a function of temperature 
for different fields as indicated. Solid line represents the calcu- 
lation described in the main text. Inset shows the comparison 
with 29 Si-NMR (I = 1/2) results on CeCu 2 Si 2 (data taken 
fromRef. 19]). (1/TiT) is normalized by {A hi ,i BO Jnmr) 2 ■ 



nated by the 4/ contribution. No significant frequency 
or field dependence of T\ could be detected in CeFePO, 
in contrast to the behavior of K(T). In Kondo lattice 
systems, \jT\T is dominated by the contribution of the 
4/ electrons, which is within simple approximations pro- 
portional to the ratio of the static susceptibility x(T) 
and the dynamical relaxation rate T(T) of the 4/ elec- 
trons, 1/TiT cx x(T)/r(T) 0, [H, Gj, [H. Further 
on, for a Kondo system one may expect T oc \pT for 
T > Tk [1(3], but merging into a constant value cx Tr- 
at low T. Since for T > 10 K, X (T) and K(T) are 
proportional, we calculated the temperature dependence 
of 1/TiT using 1/TiT = a ■ K(T)/VT, where a is just 
a T-independent parameter to scale the whole curve to 
the absolute experimental values. We added a small T- 
independent contribution (TiT) -1 = 0.12 (sK) -1 to ac- 
count for the contribution of the conduction electrons 
at high T, where the contribution of the 4/ electrons 
is vanishing. The results of this calculation shown as 
the solid line in Fig. [4j agrees very well with the data 
from 300 K down to 10 K, showing that the dynamical 
relaxation rate T indeed follows the expected y/T behav- 
ior. The deviation below 10 K is as expected, since T(T) 
should merge in a constant value instead of decreasing 
further. Further on, we compare in the inset of Fig. 0] 
31 (l/TiT) of CeFePO with 29 Si-NMR (I = 1/2) results 
of the prototypical heavy fermion metal CeCu2Si2 [H[ 
which has a similar Kondo scale. For CeCu2Si2 we specif- 
ically chose the NMR results on 29 Si in order to compare 
data of two nuclei with 7=1/2 and used the results of 
Aarts et ai, because these data were obtained at a field 
comparable to our measurements at 76.4 MHz (4.47 T). 
To account for different NMR nuclei and hyperfine cou- 
plings, (l/^T) values are divided by (A h o so • 7 nmr ) 2 . 
The comparison can then be performed on an absolute 
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scale without having introduced model dependent param- 
eters. Both curves share the same overall temperature 
dependence but CeFePO shows larger absolute values. 
The very similar T-dependences oc x(T)/vT found for 
CeFePO and CeCu2Si2 hint at similar spin-fluctuation 
relaxation mechanisms in both compounds and confirm 
that the strong correlation effects in CeFePO originate 
from the Ce-4/ electrons rather than from Fe-3d elec- 
trons. Further on, since this similarity is observed on the 
same absolute temperature scale, it confirms comparable 
Kondo temperatures in both compounds. The larger val- 
ues of the normalized \jT\T in CeFePO perfectly agree 
with the higher susceptibility. For CeFePO, x(2 K) 
saturates at 660 x 10 -9 m 3 /mol in high fields whereas 
for CeCu 2 Si 2 x(2 K) w 100 x 10~ 9 m 3 /mol was pub- 
lished L 2l| . Since both compounds have almost the same 
Sommerfeld coefficient in their electronic specific-heat, 
a larger susceptibility implies an enhanced Sommcrfcld- 
Wilson ratio R — \ll ' Ro,Ro = ^^sA^oMes) m Ce- 
FePO. With n c g = 1.73 corresponding to an average 
value Msa t = 1 MB expected for a Ce 3+ doublet ground 
state |22j, we obtain R = 5.5, which highlights the pres- 
ence of dominating FM correlations in CeFePO. This sce- 
nario is supported by the analysis of the Korringa prod- 
uct, given by T X TK 2 = S a = ^ ff /37r?ifcs7^. Magnetic 
correlations modify this relation by introducing an en- 
hancement factor K(a), given by T\TK 2 = Sa/K(a) 
where K(a) > 1 and < K{a) < 1 reflects the existence 
of AFM and FM correlations among the quasiparticles, 
respectively. For CeFePO, a value of K(a) ~ 0.065 is 
found at 2 K (4.47 T), in fact indicating strongly domi- 
nating FM correlations. A value K(a) <C 1 is also found 
in other heavy fermion metals with dominating FM cor- 
relations like YbRh 2 Si 2 Q and CeRu 2 Si 2 [24j. In con- 
trast, for CeCu 2 Si 2 with dominating AFM correlations, 
K(a) > 1 was determined [Hj|. 



To conclude, we prepared polycrystalline samples 
of CeFePO by a Sn flux technique and determined 
its ground-state properties by means of susceptibility, 
specific-heat, electrical resistivity and NMR measure- 
ments. x(T) follows a Curie- Weiss law at high tempera- 
tures with an effective moment expected for free trivalent 
Ce 3+ but levels off below 4 K at a constant, enhanced 
Pauli susceptibility xo = 660 x 10" 9 m 3 /mol. C 4f {T)/T 
increases logarithmically below 10 K and also saturates 
below 2 K at 7 = 700 mJ/K 2 mol. The temperature de- 
pendence of p(T) presents the typical behavior of a heavy 
fermion metal with a strong decrease when cooling down 
to below 30 K due to the onset of coherent Kondo scatter- 
ing. In the NMR measurements, \jT{F increases roughly 
with x(T') / VT upon cooling from 300 K to 10 K, but sat- 
urates below 10 K at a strongly enhanced Korringa term. 
Thus, CeFePO presents the typical properties of a classi- 
cal Ce-based heavy fermion metal. At low T, we observe 
the behavior expected for a Landau-Fermi liquid: a T- 
independent susceptibility xoj a T-independent C(T)/T, 
a T-independent \jT\T as well as a resistivity increment 
with T 2 . None of the investigated properties gives evi- 
dence for magnetic ordering above 0.4 K. Therefore, Ce- 
FePO may be classified as a paramagnetic heavy fermion 
metal. The absence of significant Fe-contribution to the 
effective moment at elevated temperatures indicates that 
the magnetism in CeFePO is dominated by the Ce-4/ 
electrons. Thus, the observed strong electronic correla- 
tion effects originate from the Ce-4/ electrons rather than 
the Fe-3d electrons. An enhanced Sommerfeld- Wilson ra- 
tio as well as a reduced Korringa product highlight the 
presence of dominating FM correlations. Therefore, Ce- 
FePO is likely to be located on the non-magnetic side 
close to a FM instability. 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 

J. Am. Chem. Soc. 130, 3296 (2008). 
[2] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and 

D. F. Fang, Nature 453, 761 (2008). 
[3] Z.-A. Ren, J. Yang, W. Lu, W. Yi, G.-C. Che, X.-L. 

Dong, L.-L. Sun, and Z.-X. Zhao, Europhys. Lett. 82, 

57002 (2008). 

[4] G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, 

P. Zheng, J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 

100, 247002 (2008). 
[5] C. Krellner, N. S. Kini, E. M. Briining, K. Koch, H. Ros- 

ner, M. Nicklas, M. Baenitz, and C. Geibel, Phys. Rev. 

B 76, 104418 (2007). 
[6] C. Krellner and C. Geibel, J. Cryst. Growth 310, 1875 

(2008). 

[7] C. Krellner, T. Forster, H. Jeevan, C. Geibel, and 
J. Sichelschmidt, Phys. Rev. Lett. 100, 066401 (2008). 

[8] T. M. McQueen, M. Regulacio, A. J. Williams, Q. Huang, 
J. W. Lynn, Y. S. Hor, D. V. West, M. A. Green, and 
R. J. Cava, Phys. Rev. B 78, 024521 (2008). 



[9] M. G. Kanatzidis, R. Pottgen, and W. Jeitschko, Angew. 

Chemie 117, 7156 (2005). 
[10] B. Zimmer, W. Jeitschko, J. H. Albering, R. Glaum, and 

M. Reehuis, J. Alloys Compd. 229, 238 (1995). 
[11] F. Steglich, J. Aarts, C. D. Bredl, W. Lieke, D. Meschede, 

W. Franz, and H. Schafer, Phys. Rev. Lett. 43, 1892 

(1979). 

[12] E. Dormann, in Handbook on the Physics and Chemistry 
of Rare Earths, edited by K. A. Gschneidner Jr. and 
L. Eyring (North-Holland, Amsterdam, 1991), vol. 14, 
chap. 94, pp. 63-161. 

[13] Y. Nakai, K. Ishida, Y. Kamihara, M. Hirano, and 
H. Hosono, J. Phys. Soc. Jpn 77, 073701 (2008). 

[14] H.-J. Grafe, D. Paar, G. Lang, N. J. Curro, G. Behr, 
J. Werner, J. Hamann-Borrero, C. Hess, N. Leps, R. Klin- 
geler, et al., Phys. Rev. Lett. 101, 047003 (2008). 

[15] H. Nakamura, M. Shiga, Y. Kitaoka, K. Asayama, and 
K. Yoshimura, J. Phys. Soc. Jpn 65, Suppl. B, 168 
(1996). 

[16] Y. Kuramoto and Y. Kitaoka, Dynamics of Heavy Elec- 



■5 



irons (Oxford Science Publications, 2000). 
[17] N. Biittgen, R. Bohmer, A. Krimmel, and A. Loidl, Phys. 

Rev. B 53, 5557 (1996). 
[18] D. E. MacLaughlin, Hyperfine Interactions 49, 43 (1989). 
[19] J. Aarts, F. deBoer, and D. E. MacLaughlin, Physica B 

& C 121, 162 (1983). 
[20] D. L. Cox, N. E. Bickers, and J. W. Wilkins, J. Appl. 

Phys. 57, 3166 (1985). 
[21] B. Sales and R. Viswanathan, J. Low Temp. Phys. 23, 

449 (1976). 



[22] G. Fischer and A. Herr, Phys. Status Solidi B 159, K23 
(1990). 

[23] K. Ishida, K. Okamoto, Y. Kawasaki, Y. Kitaoka, 

O. Trovarelli, C. Geibel, and F. Steglich, Phys. Rev. Lett. 

89, 107202 (2002). 
[24] Y. Kitaoka, H. Arimoto, Y. Kohori, and K. Asayama, J. 

Phys. Soc. Jpn. 54, 3236 (1985). 
[25] H. Tou, K. Ishida, and Y. Kitaoka, J. Phys. Soc. Jpn. 

74, 1245 (2005). 



